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Abstract — We consider the quantization of a transmit beam- 
forming vector in multiantenna channels and of a signature vec- 
tor in code division multiple access (CDMA) systems. Assuming 
perfect channel knowledge, the receiver selects for a transmitter 
the vector that maximizes the performance from a random vector 
quantization (RVQ) codebook, which consists of independent 
isotropically distributed unit-norm vectors. The quantized vector 
is then relayed to the transmitter via a rate-limited feedback 
channel. The RVQ codebook requires an exhaustive search to 
locate the selected entry. To reduce the search complexity, we 
apply generalized Lloyd or fc-dimensional (kd)-tree algorithms to 
organize RVQ entries into a tree. In examples shown, the search 
complexity of tree-structured (TS) RVQ can be a few orders of 
magnitude less than that of the unstructured RVQ for the same 
performance. We also derive the performance approximation for 
TS-RVQ in a large system limit, which predicts the performance 
of a moderate-size system very well. 

Index Terms — Signature quantization, tree-structured code- 
book, CDMA, MIMO, random vector quantization, generalized 
Lloyd algorithm, kd tree. 



I. Introduction 

Channel information at both the transmitter and receiver 
can increase the performance of wireless systems significantly. 
With channel information, the transmitter can adapt its trans- 
mit power and waveform to a dynamically fading channel 
while the receiver can detect transmit symbols from received 
signals. Typically, channel information can be estimated at 
the receiver from pilot signals during a training period. The 
transmitter on the other hand is usually not able to directly 
estimate a forward channel, especially in a frequency division 
duplex where a channel in one direction is not reciprocal to 
that in the opposite direction. Thus, the transmitter has to rely 
on the receiver for channel information. Normally, the receiver 
relays channel information to the transmitter via a rate-limited 
feedback channel. 

The receiver can directly quantize channel coefficients 
and feeds back the quantized coefficients to the transmitter, 
which adapts its transmission, accordingly fTl. Alternatively, 
the receiver computes and quantizes the optimal transmit 
coefficients. References 0-1161 proposed quantization of the 
transmit preceding matrix, which consists of transmit antenna 
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weights in a multiantenna channel while fT\, fS] considered 
quantization of the signature vector in CDMA. Comparing the 
two approaches, [1| showed that quantizing transmit coeffi- 
cients performs much better than direct channel quantization. 

Most of the proposed quantization codebooks require an 
exhaustive search to locate the quantized signature vector. 
The search complexity depends on the number of entries 
in the quantization codebook, which grows exponentially 
with available feedback bits denoted by B. Reducing the 
search complexity while not jeopardizing the performance is 
desirable. In [TJ, the number of coefficients to be quantized 
was reduced by projecting the signature vector onto a lower 
dimensional subspace and the coefficients were then scalar 
quantized. Although the complexity of scalar quantization is 
much less than that of vector quantization, the performance of 
the scalar quantization scheme suffers greatly. References [9|, 
[lO] proposed a search algorithm for the QAM codebook, 
which is based on a noncoherent detection algorithm ifTTl . 
The associated complexity grows linearly with B with small 
performance loss. However, the algorithm is not applicable 
when B is smaller than the number of coefficients. To reduce 
the search complexity, we propose to organize an unstructured 
random vector quantization (RVQ) codebook into a tree by 
either a generalized Lloyd (GLA) algorithm IIT2I . llT3l or 
a /-C-dimensional (kd) tree algorithm |14), (15). The tree- 
building process is computationally complex, but can be 
performed offline. An RVQ codebook was proposed in 16], Q 
and contained independent isotropically distributed unit-norm 
vectors. For moderate-size systems, RVQ performs close to the 
optimum codebook designed for a channel with independent 
identically distributed gains. In a large system limit in which 
system parameters tend to infinity with fixed ratios, RVQ 
is optimal jS), ||7] (i.e., maximizing capacity or minimizing 
interference power). 

In this work, we consider quantizing both the transmit beam- 
forming vector in multi-input multi-output (MIMO) channel 
and the signature sequence in reverse-link CDMA. With a 
tree-structured (TS) RVQ codebook, the receiver searches the 
tree for the vector that is closest in Euclidean distance to the 
optimal unquantized signature vector, which is the eigenvector 
of the received covariance matrix. We note that the obtained 
quantized vector may not give the optimal performance (i.e, 
maximizing received power in a MIMO channel or minimizing 
the interference power in CDMA). However, the performance 
loss incurred is small when B is large, while the corresponding 
search complexity is reduced by a few orders of magnitude 
from that of a full search. For smaller B, the performance, 
however, takes a significant loss. So, we modify the kd-tree 
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search algorithm to narrow the performance gap. Analyzing 
the performance of the TS-RVQ codebook is not tractable. 
Hence, we derive the performance approximation of the TS- 
RVQ codebook in the large system limit. The derived approxi- 
mation is a function of the number of feedback bits per degree 
of freedom and the normalized load and is shown to predict 
the performance of a system with moderate-to-large size very 
well. 

II. System Model 
We are interested in two wireless channel models as follows. 

A. Point-to-Point Multiantenna Channel 

We first consider a discrete-time point-to-point channel with 
Nt transmit antennas and Nr receive antennas. The A'^ x 1 
received vector is given by 



Hvh + n 



(1) 



where H = [/in^.nt] is an Nr x Nt channel matrix, v is 
an A't X 1 beamforming unit-norm vector, & is a transmitted 
symbol, n is an A'^^ x 1 additive white Gaussian noise vector 
with zero mean and co variance cr^/, and / is an identity 
matrix. Assuming an ideal rich scattering environment, we 



model h,. 



as a complex Gaussian random variable with 



zero mean and variance 1 /Nr- Thus, the received fading power 
is given by YlirJ=i ^[\^n^,nt\'^] = 1 for a given rifth transmit 
antenna. An ergodic channel capacity is given by 

C{v) = EH[\og{l + pv'^H'^Hv)] (2) 

where expectation is over distribution of H, and p = 
is the background signal-to-noise ratio (SNR). We note that 
the capacity (|2|i is a function of the beamforming vector v 
and that only a rank-one transmit beamforming is considered. 
Quantization of an arbitrary-rank transmit precoding matrix 
was studied by 16) • 

B. Reverse-Link CDMA 

We also examine a reverse-link CDMA with K users and 
processing gain N . User k is assigned the x 1 signature 
vector Sk for \ < k < K. The user signal is assumed to 
traverse L fading paths. We denote the N x N channel matrix 
for user k by 

hk,i ... 



hk,L 




ll'k,! 
hk,L 







hk.i 



hk.L 








hkA 



hk 



k,L 














hk,i 



(3) 



where fading gains for user fc, hk.i, ■ ■ ■ , hk,L are independent 
complex Gaussian random variables with zero mean and 
variance E\hk.i\'^, ■ ■■ , -B|^/c,Lp, respectively. Here we assume 
that the symbol duration is much longer than the delay spread 



and thus, inter-symbol interference is negligible. For an ideal 
nonfading channel, Ck = I for all k. 

The N X 1 received vector at the base station is given by 



where b — [bi &2 ■ 
transmitted symbols and 



r = Hb + n (4) 
. . bx] is the vector of K users' 



,Cks 



(5) 



In this work, we only consider single-user signature quan- 
tization. Without loss of generality, we assume that user I's 
signature is quantized while other signatures do not change. 
Signature quantization for multiple users was considered 
by |16|. Assuming a matched-filter receiver for user 1 



Cl 



ICiSi 



the interference power for user 1 is given by 



K 



h = y^{s\c\CkSkf = s\clHlHiCiSi 



k=2 



(6) 



(7) 



where Hi [C2S2 C3S3 . . . CkSk] is the N x (K - 1) 

interfering signature matrix. Hence, the associated output 

signal-to-interference plus noise ratio (SINR) for user 1 is 
given by 



(8) 



Similar to MIMO channels, the performance of the CDMA 
user (O also depends on the transmit signature vector Si. 

With channel information, the receiver in the MIMO chan- 
nel selects V that maximizes the capacity from an RVQ 
codebook, which is known a priori at both the transmitter 
and receiver The RVQ codebook contains independent and 
isotropically distributed unit-norm vectors and is denoted by 



(9) 



where B denotes available feedback bits. RVQ was proposed 
by ID, Q and was shown to be optimal in a large system 
hmit in which {Nt,Nr,B) — > 00. For CDMA, the receiver 
selects Si that minimizes the interference power for user 1. 

Maximizing the capacity for the MIMO channel in (O is 
equivalent to maximizing the received signal power Thus, the 
receiver selects 



V ~ are' max v'-H^Hvi 

for given H. In CDMA, the receiver selects 
Si = arg min v\CiH\HiC\v 



(10) 



(11) 



for given Hi and Ci to minimize the interference power As B 
increases, the number of entries in the codebook increases and 
so does the performance with the selected signature vector in 
both ( [Tol l and (fTTI) . However, the search complexity also grows 
since the RVQ codebook requires a full search to locate the 
optimal entry and thus, complexity can pose a serious problem 
for very large B. 
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III. Nearest-Neighbor Search 

Let SS^ denote a covariance matrix where S denotes an 
N X K effective signature matrix for both MIMO and CDMA 
channelsQ Performing a singular value decomposition (SVD) 
on the covariance matrix gives 

N 

55t = ^A»| (12) 

where Ai < A2 < • • • < Xn are the ordered eigenvalues and 
Ui is the corresponding zth eigenvector It is well known that 
the first eigenvector maximizes the quadratic form v'' SS'^v 
as follows 

SS'fy (13) 

V 

while the A^th eigenvector minimizes the quadratic form as 
follows 

Mat = argmmv^ SS^v. (14) 

V 

To quantize these optimal eigenvectors for the transmitter, the 
receiver may search the given codebook for the vector that is 
closest in Euclidean distance to either Mi or un (finding the 
nearest neighbor). To quantize Mi, the receiver selects 

Vi = arg min ||mi — Vj\\^ — argmax^{u\vj} (15) 

where is the real part of z. The right-hand side of 

(flST l follows since \\vj\\ = ||mi|| = 1. We note that Vi 
needs to be updated whenever the channel covariance SS^ 
and consequently, its eigenvectors are changed. The associated 
performance is given by 

In,,^^v\SS^Vi. (16) 

Similarly, the receiver can also quantize mat by selecting 

Vftf = argmaxSRlMjyiJj} (17) 

with associated performance 

Imin^vlSS^VN- (18) 

Finding the nearest neighbor is a classical vector quanti- 
zation problem to which there are many solutions IITtI (see 
references therein). We note that v is suboptimal and may 
not necessarily maximizes or minimizes the quadratic form 
(the received signal power in (fTOl i or the interference power 
in (fTTTl). However, the performance difference is minimal with 
a large codebook (large B). To avoid an exhaustive search to 
find the nearest neighbor, we propose to organize the RVQ 
codebook entries into a tree by applying either generalized 
Lloyd or kd-tree algorithms. 

ipor MIMO, SSt = fftff while for CDMA, 55+ = cIhIHiCi. 



A. Generalized Lloyd Algorithm 

We start with RVQ codebook V with 2^ entries. To build a 
binary tree, we iteratively divide the RVQ entries into two 
groups with the generalized Lloyd algorithm (GLA) [17|. 
We note that with this algorithm, the tree produced may not 
be balanced. Building the tree becomes more complex as B 
increases. However, it should not incur any additional delay 
since the tree can be produced offline. 

To find the optimal vector v, we apply the encoding method 
proposed by ifTSl . Since the method is only applicable to a 
real vector, we transform an x 1 complex vector Vj into 
a 2N X 1 real vector [vj^, I'JJ"^, where Vj^r = ^{vj}, 
Vj.i — ^{vj}, and is the imaginary part of z. The 

stated encoding method produces v that is closest in Euclidean 
distance to the desired eigenvector with search complexity 
growing approximately lineariy with B. 

B. Kd-Tree Algorithm 

A kd-tree is a data structure for storing points in a k- 
dimensional space. The kd-tree algorithm produces an unbal- 
anced binary-search tree by clustering the codebook entries by 
dimension at each step [14J, yj]. Since the kd-tree algorithm 
searches for the nearest neighbor of a real vector, we again 
transform an iV x 1 complex eigenvector into a 2N x 1 
real vector before quantization. First, we construct an RVQ 
codebook with 2N x 1 real unit-norm vectors and store the 
codebook at the root node. Then, we calculate a median for 
the first elements of all vectors in the codebook and select the 
element closest to the median as the pivot. By comparing the 
first element of the vector to the pivot, all vectors are divided 
into 2 groups, which will be stored in the left-and right-child 
nodes. Next we move to either the left- or right-child nodes 
and compute the pivot for the second dimension to divide the 
vectors into 2 groups for that node. We operate on the next 
dimension each time we move down the tree and iterate the 
process until each node contains only one vector. The detailed 
steps for building kd-tree are shown in llT4l . llTSll . Building 
the kd-tree is relatively faster than building the tree by GLA 
since we examine one dimension at a time and this can be 
done offline. 

To quantize the desired vector, we start at the root node. 
By comparing the element of the vector with the pivot of 
the present node in the specified dimension, we move down 
to either the left- or right-child nodes. We transverse the 
tree until the leaf node is reached and the candidate entry is 
obtained. The algorithm proposed in [141, [151 makes certain 
that the candidate is the nearest neighbor by comparing the 
distance of the candidate and the vector to quantize with 
that of other nodes. The kd-tree search produces v that is 
closest in Euclidean distance to the desired eigenvector among 
the entries in the RVQ codebook. The associated search 
complexity is a fraction of that of a full search as simulation 
results wiU demonstrate. 

C. Performance Approximation 

Evaluating the quadratic form with the nearest-neighbor 
quantized vector is an open and difficult problem. Hence, we 
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approximate the nearest-neighbor criterion with the close st-in- 
angle one. For a closest-in-angle search, the receiver selects 

Vi = argmax ImId^P ~ argmax cos^ (pj (19) 

where (j>j is the angle between Ui and Vj, and the correspond- 
ing performance is given by 



/max = vlSS^'Vi 



N 



t...|2 



Xi\v{u, 

i=2 



(20) 



(21) 



where SVD in is applied. (Similarly to ([19) and (EO), we 
can also define vn and /min-) 

Evaluating the expectation for /,nax for finite N, K, and 
B is not tractable. Hence we analyze /max in a large system 
limit where [N, K,B) ^ oo with fixed ratios. The first term 
on the right-hand side in ( 1211 1 was shown by |6| to converge 
as follows 

\vlui\^ ^1-2-^ (22) 

almost surely as {N,B) ^ oo with fixed B = B/N. With 
infinite feedback {B — oo), Vi — > Ui and hence, I^ImiP — > 1. 
On the other hand, with no feedback, Vi is randomly selected 
and Itjlitip — > 0. The second term in ( l2Tl i converges to the 
following limit. 

Lemma 1: As {N, B, K) ^ oo with fixed B = B/N and 
K — K/N, we have 



N 
i=2 



^9ssi (A) dA 



(23) 



assuming that the eigenvalue density for SS'^ converges to a 
deterministic function gss^i )- 

Proof: Since entries in the RVQ codebook are uniformly 
distributed on the unit-norm hypersphere, it was shown by f4l 
that Prllv^Mip > s^} is proportional to the surface area of 
the spherical cap of the A^-dimensional hypersphere described 
by 



N 

El 

n=l 



(24) 



intersecting with > s\, where w„ is the n\h complex 

entry of Ui and r denotes the radius0 The surface area of the 
described spherical cap is given by ||4l 



AN{r,si) 



2ti 



N 



{N-l) 



2\N-1 



(25) 



We would like to evaluate Pr{|'i)^^M2p > conditioned 
that {v''ui)^ > sf. Similar to the results in ||4l, we can deduce 
that the conditional probability is proportional to the surface 
area of the intersection between the two spherical caps. The 
volume of the intersection is given by 



Ml e 



N 
n=l 



r, rj > sj and > s^j 
(26) 



Later we will set r = 1. 



and can be computed with spherical coordinates as follows 



VN^r, 81,82) = (27r) 



N 



ri=si Jr2=S2 



rs • • • r„ dra • • • drN rir2 dr2dri . (27) 



We note that the multiple integral in the brackets in (|27T i is the 
volume of an {N — 2) -dimensional hypersphere with scaling 
factor of (27r)^-2 g]. Thus, 



VNir, 81,82) 



(27r) 



N 



2^-2(7V-2)7,^^,^ 
X (r-r^-r2)^-2^r2dri 
{2n)^ 



rir2 



{r~s\-8l)^ 



2^N 



(28) 
(29) 



The associated surface area is obtained by differentiating the 
volume VN{r, s^, s"^) with respect to r and is given by 

<WN[r, 81,82) 



AN{r,8i,82) 



dr 



27r 



N 



-Ar'-sl-8l)^'' 



(N-l) 

Thus, the conditional probability 

PT{\v^U2f>sl I \i^Ui\^>sl} 

_ Pr{|'MtM2p > 8^, |i)tM^|2 > sj} 
~ Pr{|i}tMi|2 > sj} 

_ ANir,8l,S2) 



(30) 
(31) 



AN{r,8i) 



1 



p2 — c2 



W-1 



(32) 
(33) 

(34) 



and the corresponding cumulative distribution function (cdf) 
for |t}|M2p given that [■ujitip > and r = 1 is given by 

N-l 

l-r,t„. |2 ^ ,2^l _ 1 _ ^ 1 ^ \ _ ^^2,5) 



For the ith eigenvector where 17^!, the cdf for Itjlwip is the 
same as shown in (1351 ). Thus, with the cdf, we can compute 
the conditional expectation as follows 



E[\vlu.\^\\vlu,\'>8i]^-{l-8i) (36) 

which is converging to zero as — > 00. Pr{|{j|u.ip > s^} 
depends on the number of entries in the RVQ codebook or 
available feedback bits. As {N, B) ^ 00, fF\ has shown that 



1 



1 - sf ^ 2- 



(37) 

If the eigenvalue density of SS^ converges to a deterministic 
function ffsst(-). 



1 ^ r 



Agsst (A) dA. 



(38) 



as {N, K) 00 with fixed K/N. Combining ^ and ( [38] |, 
we have Lemma [T] ■ 
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From Lemma [T] as B increases, Vi and Ui are becoming 
perpendicular and |t3|i6ip 0. We apply Lemma[T]to obtain 
the following performance approximations for the nearest- 
neighbor search for MIMO and CDMA models in section HIl 

Theorem 1: As {Nr,Nt,B) oo with fixed B/Nt and 
Nr = Nr/Nt, the capacity of MIMO channel with the nearest- 
neighbor beamforming vector can be approximated as follows. 




(1-2- 



(39) 



(40) 



Proof: We evaluate the received power vIh'^Hvi by 
applying Lemma [T] and (|22] |. Thus, 



XgninWdX (41) 



H^H is a Wishart matrix with the well known asymptotic 
eigenvalue distribution. In the large system limit, the maximum 
eigenvalue converges to |,19J 



maxV 



H^H) 



1 



'Nr 



and 



(42) 



(43) 



where Qh^hW is given by |19|. Combining the capacity 
expression in (|2]i and (|4T1 |-(|43]| gives Theorem [T| ■ 
Theorem 2: For CDMA, suppose the channel gain for user 



-¥ a while channel gains for interfering 

2 ^ 1 for 2 < A: < X as (TV, K) oo. 



1. T.r=iE\hu? ' 
users, X^fe^i E\hkA 
The large system SINR for user 1 with the nearest-neighbor 
quantized signature and a matched filter is approximated as 
follows. 



7(t;) « 7(i}) 



(44) 



(1-^)2(1-2- 



■.Q<K <1 
: K >1 



(45) 



The proof is similar to that for Theorem [T] with the asymptotic 
minimum eigenvalue and eigenvalue density for the inter- 
ference covariance CiHIHiCI shown in |[T9l, HqI. If the 
channel gains across interfering users are not uniform (i.e., 
"^^^i E\hk^i\^ is arbitrary for fc), we apply the asymptotic 
eigenvalue density of CiHIHiCI with a non-uniform power 
allocation from 1191 . ||20 |. In section |V] we will compare 
the asymptotic approximations derived here with simulation 
results. 



IV. Modified Kd-Tree Search 

We modify the search for the kd-tree so that we transverse in 
the direction that maximizes or minimizes the quadratic form 
v'^ SS^v. Suppose v^r) = ^{v} and — ^{v}. Expanding 
the quadratic form gives 

+ vl^^{SS^v^r) + vl^^iSS^v^,-, (46) 

where we use the fact that SS^ is Hermitian. We note that 
the two cross terms in (l4&i are much smaller than the first 
two quadratic terms. Hence, our objective is to maximize 
or minimize the two quadratic terms in (l46t . which only 
depend on JRjS'S't} . To quantize t^j-^j, we start with the RVQ 
codebook whose entry has norm 1/2. Similar to the nearest- 
neighbor search, we start at the first dimension and determine 
the pivot. Then, we divide all entries into 2 child nodes and 
compute the second pivot for each node. We move to the child 
node whose pivot gives the larger value of t)^^.^ SKIS' S'^}t>(r)- 
This differs from the nearest-neighbor search, which does not 
take SRjS'S'^} into account. We continue the process until we 
reach the leaf node and hence, produce the candidate vector 
The steps are summarized in Algorithm [T] and Algorithm |2] 
which compares the candidate with the surrounding nodes. 

Algorithm 1 Modified kd-tree search 
1: Start at the root of a kd-tree codebook with the channel 

covariance matrix SS^ . 
2: Set n — 1 where n is the index of element in an A^- 
dimensional vector, 
while Not a leaf node do 
if n > TV then 
n = l 



end if 

pi"' = [0 . . . PL.n+l . 

vector of the left-child node 

= [0 . . . PR.n+l . 

vector of the right-child node. 

>^ > p^^^^SS^p^^^ then 
Move to the left-child node. 











0]t where is the 
0]^ where pn is the 



if Pi ss^pI 



else 

Move to the right-child node, 
end if 
n = n + I 
end while 

V = vector from the leaf node. 

n = 1 

s = compare(root node, n, S, I, v) 
return s 



This modified kd-tree search does not produce the optimal 
vector that maximizes or minimizes the quadratic form, but 
can perform close to the optimum when N is large as shown 
in numerical examples. Analyzing performance of this scheme 
is not tractable. However, the performance is upper bounded 
by that of the full search, which was analyzed in a large system 
limit by 0, lED. 
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Algorithm 2 function: s = compare(node, n, S, I, v) 

1: if Leaf node, then 

2: Sc — vector from the leaf node. 

3: Ic = slSSU, 

4: it Ic> I then 

5: I ^ Ic 

6: S = Sc 

1: return s 

8: end if 

9: end if 

10: pi"'' = [0 . . . Pc,n • . • 0]^ where Pc is the vector 
of the current node. 

11: /p=pi")^SStp(") 

12: if /p > / then 

13: if n > then 

14: n = 1 

15: end if 

16: if p^fsS^p^i^^ < p%'^ SS^p%'^ then 

17: n = n + 1 

18: compare(left-child node, n, 5', /, s) 

19: else 

20: n = n + 1 

21: compare(right-child node, n, S, I, s) 

22: end if 

23: else 

24: n = n + 1 

25: compare(left-child node, n, /, s) 

26: compare(right-child node, n, S, I, s) 

27: end if 



V. Numerical Results 

Fig. [T] shows the channel capacity of the quantized beam- 
former for the MIMO channel with TV* = 3 and iV^ = 4 and 
p = 10 dB. We compare different quantization schemes and 
note that RVQ with full search gives the maximum capacity 
for a given feedback. When there is no feedback (B — 0), 
the transmitter deploys a random beamformer and thus, all 
schemes give the same performance. As available feedback in- 
creases, the performance of the quantized beamforming vector 
increases with different rates depending on which quantization 
scheme is used. We remark that the proposed kd-tree with a 
modified search performs close to RVQ for small B. We also 
plot the performance of the GLA and kd-tree schemes from 
Section |III] and also that of the suboptimal QAM codebook 
proposed by [9|. Reference [9| also proposed the optimal 
QAM codebook scheme that searches for the codebook entry 
closest in angle to the optimal eigenvector Here we opt to 
compare our proposed schemes with the suboptimal QAM 
codebook since the suboptimal QAM codebook performs very 
close to the optimal QAM codebook with significantly less 
complexity. In Fig. [T] we see that quantization of the eigen- 
vector does not perform well for small B, but its performance 
is very close to that of RVQ for large B. The QAM codebook 
performs a bit better than the other two schemes, but is not 
appHcable for B < 2Nt. 

In addition to channel capacity, we also examine the com- 



N =3;N =4;SNR = 10dB 

t ' r ' 




Number of feedback bits 

Fig. 1. Shown is the ergodic capacity of the quantized beamforming vector 
for MIMO channel with 3 transmit and 4 receive antennas. The performance 
of different quantization schemes are shown with B feedback bits. 

putational complexity of each quantization scheme required to 
search for the selected entry in the codebook. Fig. |2]shows the 
average number of matrix inner product computations used for 
each scheme versus the capacity with the same set of param- 
eters as the previous figure. As expected, the RVQ codebook 
with a full search requires the most number of inner product 
computations, which grows exponentially as B increases. At 
a capacity of 4.2 bits per channel use, the kd-tree requires 
almost two orders of magnitude less computations than RVQ 
does while the kd-tree with the modified search requires almost 
one order of magnitude less. The QAM codebook is the least 
complex for a large capacity, however it is not available for a 
low capacity. 



N =3;N =4;SNR = 10dB 

t r 




3.4 3.6 3.8 4 4.2 4.4 



Capacity (bits/channel use) 

Fig. 2. The average number of matrix inner product computations required to 
locate the selected entry for various quantization codebooks is shown with the 
ergodic capacity for MIMO channel with Nt = 3, A^r = 4, and p = 10 dB. 

We compare the large system capacity approximation in 
Theorem [T] with the capacity of a finite-size MIMO channel 
with Nj./Nt — 1 and show the comparison in Fig. [3] We 
note that there is a substantial gap between the theoretical 
approximation and the simulation results especially for the 
2x2 system. The gap narrows as the system size increases to 
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4x4 and is expected to narrow further as the size increases. For 
the 2x2 channel, one feedback bit per transmit antenna does 
not improve the capacity by much. However, as the system 
size increases, one feedback bit per antenna could potentially 
increase the capacity as much as 50% from the zero-feedback 
capacity. 



N/Nj = 1 ; SNR = 5dB 





/ e < 








Large sys. approx in Thm. 1 

KD-TREE w/ = 2 

KD-TREE w/N^ = 4 



B/N| 

Fig. 3. The asymptotic capacity approximation is sliown with the capacity 
of 2 X 2 and 4x4 channels. 

The simulation results for CDMA are shown in Fig. |4] The 
kd-tree codebook with a modified search has performance that 
is close to that of the RVQ codebook while the eigenvector- 
quantized schemes perform much worse, especially for small 
B. We also plot the large system approximation in Theorem|2] 
Unlike the MIMO results, the derived approximation predicts 
the performance of the finite-size CDMA quite well since the 
system size here {N — 10) is much larger than that in MIMO. 
We note that results for the RVQ codebook with large B are 
missing due to the very large memory requirement to store the 
codebook and the large computing power required to locate 
the codebook entry. In the case of the GLA and the kd-tree 
codebooks, the results for large B are also missing because 
building a tree with large number of entries also demands very 
high computing power 

We also show the large system performance of the RVQ 
codebook derived in |21|. We remark that with one feedback 
bit per processing gain, SINR for the RVQ codebook is about 
4 dB higher than that for the eigenvector-quantized codebooks. 
However, this advantage comes with a much larger search 
complexity as we will see in Fig. |5] 

Besides capacity, we would like to compare the complexity 
of each scheme to locate the selected entry from the RVQ 
codebook. The computational complexity will be measured 
by the number of inner products between two iV-dimensional 
vectors. For an exhaustive search, the number of inner-product 
computations increases exponentially with B. For the proposed 
binary tree with the nearest-neighbor search, the number of 
inner-product computations depends on the tree's depth and 
only increases linearly with B. The complexity of the kd- 
tree search is even less since in each search step, only one 
dimension of an A^-dimensional codebook entry is required. 
Thus, the search complexity of the kd-tree is proportional to 



N = 10; K = 5; L = 1; SNR = 10dB 

10 I ^ ^ ^ _ " 




0.5 1 1.5 2 2.5 3 3.5 4 

B/N 

Fig. 4. SINR with quantized signatures from different quantization schemes 
are compared. 

BjN . This also applies to the modified kd-tree search as well. 

We compare the search complexity associated with different 
codebooks whose SINR performance is shown in Fig. |5] 
Similar to the results for the MIMO channel, the average 
number of inner product computations used to search the 
selected entry is shown. The RVQ codebook is the most 
complex while the kd-tree and QAM codebooks are the least 
complex. At 9 dB, the kd-tree codebook with a modified search 
requires about 3 orders of magnitude in search complexity 
less than RVQ does. The QAM codebook also requires few 
computations, however, it is not valid for small number of 
feedback bits. 




23456789 10 
SINR (dB) 



Fig. 5. The average number of inner product computations to search for 
selected entry in the codebook for each quantization scheme is shown with 
SINR in dB. 



VI. Conclusions 

We have proposed tree-structured RVQ codebooks to quan- 
tize the transmitted signature for CDMA and the beamforming 
vector for the MIMO channel. The tree structure is gener- 
ated by either the GLA or kd-tree algorithms. For small- 
to-moderate feedback, the kd-tree codebook with a modified 
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search has performance that is close to that of the RVQ with 
the exhaustive search. Thus, our proposed scheme comple- 
ments the QAM codebook proposed by Ryan et al |9|, which 
does not perform well in a low-feedback regime. For large 
feedback, quantizing the eigenvector with the GLA or the 
kd-tree codebooks produces a performance that is close to 
RVQ. The search complexity of the proposed schemes are a 
few orders of magnitude less than that of RVQ for a given 
performance. We also approximated the performance with 
large system limit, and numerical examples have shown that 
the approximations can predict the simulation results well for 
a moderate-size system. 

In this study, we only consider quantization of a rank-one 
transmit beamforming vector in MIMO channels and single- 
user signature quantization in CDMA. To extend the existing 
results to MIMO channels with an arbitrary-rank precoding 
matrix or uplink CDMA with multiple-signature quantization, 
we can vectorize the precoding matrix or multiple signatures 
and apply the GLA or kd-tree algorithms to construct a tree. 
However, the mentioned scheme can be highly suboptimal. 
Thus, improving the scheme and analyzing the associated 
performance are interesting open problems. 
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